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We focus on a linear chain of N first-neighbor-coupled logistic maps at their edge of chaos in 
the presence of a common noise. This model, characterised by the coupling strength e and the 
noise width Omax, was recently introduced by Pluchino et al [Phys. Rev. E 87, 022910 (2013)]. 

They detected, for the time averaged returns with characteristic return time r, possible connections 
with q-Gaussians, the distributions which optimise, under appropriate constraints, the nonadditive 
entropy Sq, basis of nonextensive statistics mechanics. We have here a closer look on this model, and 
numerically obtain probability distributions which exhibit a slight asymmetry for some parameter 
values, in variance with simple g-Gaussians. Nevertheless, along many decades, the fitting with 
g-Gaussians turns out to be numerically very satisfactory for wide regions of the parameter values, 
and we illustrate how the index g evolves with {N,T,e,amax). It is nevertheless instructive on 
how careful one must be in such numerical analysis. The overall work shows that physical and/or 
biological systems that are correctly mimicked by the Pluchino et al model are thermostatistically 
related to nonextensive statistical mechanics when time-averaged relevant quantities are studied. 

PACS numbers: 05.45.Ra, 74.40.De, 87.19.1m 


I. INTRODUCTION 

Synchronization has long been observed in complex 
systems and extensively studied in the literature 0,0. In 
these studies, coupled maps are considered as an impor¬ 
tant theoretical model for these systems Q. Moreover, 
since many biological systems evolve in noisy environ¬ 
ments, it is important to analyse the effect of noise in 
such coupled maps. Recently, the effect of weak noise 
on globally couMed chaotic units has been studied for 
several systems [3,0] . 

The model that we focus on here is a linear chain of 
N coupled logistic maps, with periodic boundary condi¬ 
tions, and can be written as 


where p G [0, 2] is the map parameter and this function 
is taken in module 1 with sign in order to fold the iter¬ 
ates of the maps back into the map interval [—1,1] if the 
noise takes them out of this interval. All elements of the 
system will be kept at the chaos threshold by fixing the 
/i parameter at pc = 1.40115518909.... 

Very recently, the probability distribution functions 
(PDFs) of the returns for this model have been inves¬ 
tigated numerically by Pluchino et al.0] and fat-tailed 
distributions have been reported, which can be fitted by 
g-Gaussians. The physical quantity under investigation, 
the returns, is a commonly used one in areas such as 
turbulence 0] , finance 0 H , DNA sequences 0 , and 
earthquake dynamics |lll - ll4| in the literature. It is de¬ 
fined as 


xl^^ = {l-e)f {xl) + ^[f {xl ^) + f {xl+^)]+a{t) (1) 

where e G [0,1] is the local coupling, and a{t) is an ad¬ 
ditive random noise uniformly distributed in [0, Umaa:]) 
which fluctuates in time but is equal for all maps. Here, 
the ith logistic map at any time t is given as 

f{xl)=l-pi\xl\' ( 2 ) 
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Adt — dt+T — dt (3) 

where 

1 ^ 

= ■ ( 4 ) 

i^l 

The main result of 0 has been given in Fig. 3 of 
that paper where the PDFs of the normalized returns 
(normalized to the standard deviation of the overall se¬ 
quence) have been plotted for various amax values. We 
have reproduced here in our Fig. [1] two representative 
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normalized returns 



FIG. 1: Probability distribution functions of the normalized 
returns for N — 100 logistic maps at the edge of chaos with 
r = 32. The parameter values used here are exactly those of 
Fig. 3a (case (a)) and of Fig. 3d (case (b)) of Q. By expr = 1 
we mean that only one experiment has been done in each of 
these illustrations. 


cases, namely, amax = 0.002 and (Jmax = 0.3 to be com¬ 
pared with Fig. 3a and Fig. 3d of respectively. It 
is evident that the g-Gaussian curves with predicted q 
values are very reasonable for these examples {q = 1.54 
for (Tmoa: = 0.002 and q = 1 for cFmax = 0.3). However, 
the distributions span less than 3 decades since the time 
series that have been used are relatively short {t = 10®). 
Therefore it is interesting to check the behaviour with 
longer time series in order to span more decades of the 
PDF. This is what we have done: see Fig. [2l which has 
been obtained by considering a large number of exper¬ 
iments (250 and 500). This means that the time series 
analysed are as long as 250 x 10® and 500 x 10® respec¬ 
tively. As seen clearly from Fig. [21 when one more decade 
is exhibited for the tails, the departure from the a priori 
admissible g-Gaussians becomes clearly visible. 

In addition to this, it is also interesting to check 
whether these typical numerical PDFs are asymptotic 
in the sense of the N ^ oo limit. If this is the case. 



normalized returns 



normali z ed returns 


FIG. 2: Probability distribution functions of the normalized 
returns for N = 100 logistic maps at the edge of chaos with 
T = 32 for Omax = 0.002 (a) and Omax = 0.3 (b). In order to 
analyse the tails of the distribution, large number of experi¬ 
ments is used. The size of the time series therefore equals to 
the number of experiments times the time steps used in each 
experiment (i.e., exprxlO®). 


one would expect the PDF curves to remain basically in¬ 
variant as the size N of the system increases. We have 
checked it: see two illustrative cases in Fig. |3l It is eas¬ 
ily seen that the N = 100 size cannot be considered as 
sufficient for the analysis of the asymptotic behaviour, 
especially for small values of Umax- 

It is therefore clear at this stage that, if we are par¬ 
ticularly interested in studying the asymptotic behaviour 
of PDFs of the normalized returns of this system, both 
the size N of the system and the length of the time se¬ 
ries must be quite large. To further refine numerically 
the present analysis, we have also studied the same two 
cases with large enough value for N {N = 500000 for 
CTmax = 0.002. and N = 200000 for cTmax = 0.3), and 
with a time series of 100 x 10® steps. These conditions 
are sufficient to observe the behaviour up to almost 6 
decades. The results are given in Fig. U) The strong 
departure from the g-Gaussians with early predicted q 
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normalized returns 



normalized returns 


h = 0) as illustrated in Fig. The representation of 
the same data using the g-logarithmic function, defined 
aslngy= (yi“«-l)/(l-g) {y > 0; Ini y = Iny), has also 
been given in Fig. [Hd. We notice that, in the asymptotic 
regime, the PDF for large values amax approaches a q- 
Gaussian with q < 1, rather than a Gaussian, as first 
advanced in Ref.@ far from the asymptotic regime: see 
Figs. 01: and 011. 


II. 5-GAUSSIAN approximants 

Let us now attempt a wider look at the same prob¬ 
lem, by extending the preliminary results given in [^. 
A simple y-Gaussian approximant describes fairly well 
the numerical data unless we focus on a parameter value 
where an asymmetry emerges in the PDF. The index q 
is generically expected to depend on (e,CTmax)T). 

First of all, we need to define the procedure that we 
are employing in order to predict the best value of q pa¬ 
rameter of the related y-Gaussian at each tuple. The 
y-Gaussians can be obtained by optimising the entropy 
Sq = ^1 - /(? “ 1) [Hi and can be given for 

—oo < g < 3 as 

P = Aq^qe-^^-" (7) 

where e^, known as ^-exponential, is defined by 


FIG. 3: Probability distribution functions of the normalized 
returns for N logistic maps at the edge of chaos with r = 32 
for N = 100, N = 100000 and N = 500000. It is evident that, 
as N increases, the asymptotic PDFs are attained slowly for 
small Umax values (i.e., case (a)), whereas it is quite quick for 
larger amax values (i.e., case (b)). 


gy _ / [1 +(1, i + (i-9)y>0 

0 , else. 

By inverting this function we obtain the g-logarithm 
function defined above. 

The y-dependent coefficient Aq in Eq. 0 is given by 


values becomes now very obvious. Moreover, a slight but 
neat asymmetry emerges for small amax- In order to fit 
this kind of behavior, one can use a non-symmetric ex¬ 
tension of y-Gaussian, namely given by 






q < 1 
q=l 

1 < g < 3. 


(9) 


p(x [l-{1 - q){bx‘^ + (5) 

or even 

p oc [l — (1 — q)[bx‘^ -I- gx^ -|- hx'^)] . (6) 

These expressions recover the symmetric g-Gaussian for 
g = h = 0. Strictly speaking, the case g ^ 0 and h = 0 
is mathematically inadmissible since it would lead to a 
runaway to infinity. In contrast, the case g ^ 0 and h > 0 
is perfectly admissible. However, in practice, unless we 
explore very large values of x, we can consider h = 0 
with no sensible damage. A very good fitting has been 
obtained by using q = 1.153, b = 0.655, g = 0.006 (with 


The quantity Bq characterizes the PDF width w of the 
distribution as follows: 


where w is, for q < 5/3, related to the standard deviation 
s through 

(5 — 3g)s^ = (3 — q)w‘^ . (11) 

In our simulations, there are only two free parameters, 
namely q and i?q, to be adjusted under the assumption 
that the number of decades that we are observing still 
shows a symmetric PDF. For a given set (e,crmax,T), 





























FIG. 4: Asymptotic PDFs of the normalized returns with large N and large number of time steps for amax ~ 0.002 (a) and 
<7max = 0.3 (c). For both cases r = 32. The (/-logarithms versus the square of normalized returns have already been given for 
both cases in (b) and (d) respectively. 


we adjust the best value of q by looking at the curves 
Inq(PDF) versus (normalized returns)^. The procedure 
consists in determining the value of q which provides the 
best approach to a straight line. More precisely, we fit the 
curves resulting from various values of q with a quadratic 
function 

y = a + Px + , ( 12 ) 

and choose the value of q which provides the minimal 
value for the parameter 7 . Then, for this q value, we 
calculate Aq and Bq parameters from the time series. Fi¬ 
nally, using all these obtained parameter values, we plot 
the best (/-Gaussian for the chosen set {e,<7max,T). Two 
typical examples are given in Fig.[5l For these examples, 
the minimal value of 7 , denoted by 7 mm > have been ob¬ 
tained for q = 0.964 and q = 1.022, respectively. Then 
these values of q have been used to determine the values 
of the Aq and Bq parameters. Finally, the PDFs of the 
normalized returns are approached as shown in Fig. Or 


and 5d. It is evident from the figure that the best q- 
Gaussian approximants corroborate the behaviour of the 
system for more than 6 decades in each case. The ob¬ 
tained results are summarized in Table I and plotted in 
Fig. [ 6 ] for typical values of {e.Umax)- Unless the Umax 
values are very small, the tendency is always towards a 
(/-Gaussian with (/ < 1. For very small Umax, 9 > 1 are 
also observed but for such cases the asymmetry becomes 
important. 

In our simulations, we prefer to set r = 128 since it is 
evident from Fig. [3 that if r is not very small the PDF 
remains the same as r is increased. 

We have also checked the effect of e for a fixed value 
of Umax- A clear tendency of increasing q values as e 
is decreasing has been observed with a saturation value 
around q ~ 0.965 as given in Table II and plotted in 
Fig. [H 

Finally, we have simulated the case /r = 2 without 
noise. As expected, in this case, a clear Gaussian with 6 
decades can be seen in Fig. [S] since all individual logistic 
























































FIG. 5: Two typical examples of choosing the most appropriate q value using the quadratic function for a case g < 1 (a) and for 
a case g > 1 (c). The appropriate g value is taken as the one that gives the smallest 7 value (denoted as 7 min). The probability 
density of the same examples using the g value obtained from the explained procedure for the case g < 1 (b) and for the case 
g > 1 (d). 
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FIG. 6 : The behaviour of the index g as a function of a max 
for two typical values of e. 



normalized returns 


FIG. 7: The effect of r on the behavior of the PDF. 
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e = 0.000001 

(^max 

Bq 

Bq 

g 

0.001 

0.4023 

0.5171 

1.022 

0.05 

0.3957 

0.4840 

0.978 

0.1 

0.3989 

0.4897 

0.986 

0.2 

0.3944 

0.4771 

0.968 

0.3 

0.3944 

0.4771 

0.968 

0.4 

0.3953 

0.4819 

0.975 

0.5 

0.3938 

0.4744 

0.964 

0.6 

0.3937 

0.4737 

0.963 

0.7 

0.3930 

0.4704 

0.958 

00 

O 

II 

^max 

^q\/Bq 

Bq 

g 

0.0015 

0.4214 

0.6154 

1.125 

0.02 

0.4077 

0.5441 

1.054 

0.15 

0.3912 

0.4613 

0.944 

0.2 

0.3888 

0.4494 

0.925 

0.3 

0.3926 

0.4684 

0.955 

0.4 

0.3912 

0.4613 

0.944 

0.45 

0.3885 

0.4482 

0.923 

0.5 

0.3873 

0.4423 

0.913 

0.55 

0.3950 

0.4805 

0.973 


TABLE I: For r = 128 and two typical values of e, the val¬ 
ues of the parameters Aq^jBq, Bq and q as amax varies are 
indicated. 


(^max — 0.5 

e 

^q\/Bq 

Bq 

g 

0.000001 

0.3938 

0.4744 

0.964 

0.00001 

0.3937 

0.4737 

0.963 

0.0001 

0.3941 

0.4757 

0.966 

0.001 

0.3939 

0.4751 

0.965 

0.01 

0.3939 

0.4751 

0.965 

0.1 

0.3914 

0.4625 

0.946 

0.2 

0.3895 

0.4531 

0.931 

0.5 

0.3878 

0.4446 

0.917 

0.8 

0.3873 

0.4423 

0.913 


TABLE II: For r = 128 and a typical value of amax, the values 
of the parameters Aq\jBq, Bq and q as e varies are indicated. 


maps are at their mostly chaotic points. Whenever a 
nonzero noise term is included, the system is not able to 
achieve the Gaussian due to the persistence of the noise 
on each map. 


III. CONCLUSIONS 

We have numerically studied, along many decades, the 
time-averaged returns in the model recently introduced 
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FIG. 8: The behaviour of the index q as a function of e for a 
typical value of amax- 



FIG. 9: The behaviour of the PDF for a typical case of p = 2. 


in @ by varying the number N of logistic maps in the 
chain, the first-neigbour-coupling constant e, the width 
f^max of the common noise, and the characteristic return 
time T. 

Our results can be summarised as follows. When the 
control parameter p of the logistic maps is taken to be 
close to the edge of chaos {fj-c), where the Lyapunov expo¬ 
nent vanishes, slightly asymmetric distributions are ob¬ 
served for some parameter values. Unless the number of 
decades is quite large, this asymmetry can be neglected 
and the distributions satisfactorily admit g-Gaussian ap- 
proximants even for these parameter values. We have 
illustrated the effect on the index q of the various pa¬ 
rameters of the model. 

When we choose for the parameter p of the logistic 
maps not its value at the edge of chaos but any other 
value such that the Lyapunov exponent is positive, we 
obtain numerical PDFs that, as expected, are well htted 
by simple Gaussians along six PDF decades when the 
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system has no noise. In the presense of noise, however, 
the system is unable to achieve the Gaussian PDF due 
to the persistent perturbation of the noise. 

The overall scenario is that, when we consider 
time-averaged relevant quantities such as the returns, 
Boltzmann-Gibbs statistics (with its Gaussians, consis¬ 
tent with the classical Central Limit Theorem) emerges 
when strong chaos (positive Lyapunov exponent) is 
present, and nonextensive statistical mechanics (with its 
g-Gaussians, consistent with the g-generalized Central 
Limit Theorem 0 ) emerges, or nearly emerges, when 
weak chaos (zero Lyapunov exponent) is present. These 
facts reinforce the results recently obtained with time- 
averaged quantities for the (conservative) standard map 


0 . Biological and other systems ^pear to be well mim¬ 
icked by the Pluchino et al model @ . The present results 
might be useful for examining such systems. 
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